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Abstract The precise knowledge of the atmospheric neu- 
trino fluxes is a key ingredient in the interpretation of the 
results from any atmospheric neutrino experiment. In 
the standard atmospheric neutrino data analysis, these 
fluxes are theoretical inputs obtained from sophisticated 
numerical calculations. In this work we present an alter- 
native approach to the determination of the atmospheric 
neutrino fluxes based on the direct extraction from the 
experimental data on neutrino event rates. 
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1 Introduction 

One of the most important breakthroughs in particle 
physics, and the only solid evidence for physics beyond 
the Standard Model, is the discovery - following a va- 
riety of independent experiments jllj - that neutrinos 
are massive and consequently can oscillate among their 
different flavor eigenstates. The flavor oscillation hypoth- 
esis has been supported by an impressive wealth of ex- 
perimental data, one of the most important pieces of 
evidence coming from atmospheric neutrinos. 
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Atmospheric neutrinos are originated by the colli- 
sions of cosmic rays with air nuclei in the Earth atmo- 
sphere. In the collision mostly pions (and some kaons) 
are produced and they subsequently decay into electron 
and muon neutrinos and anti-neutrinos. These neutrinos 
are observed in underground experiments using differ- 
ent techniques 1,- In particular in the last ten years, 
high precision and large statistics data has been avail- 
able from the SuperKamiokande experiment 1 which 
has clearly established the existence of a deficit in the fi- 
like atmospheric events with the expected distance and 



energy dependence from 



oscillations with oscil- 



lation parameters Am ~ 2 x 10 eV and tan 6 = 1. 
This evidence has also been confirmed by other atmo- 
spheric experiments such as MACRO and Soudan 2. 

The expected number of atmospheric neutrino events 
depends on a variety of components: the atmospheric 
neutrino fluxes, the neutrino oscillation parameters and 
the neutrino-nucleus interaction cross section. Since the 
main focus of atmospheric neutrino data interpretation 
has been the determination of neutrino oscillation pa- 
rameters, in the standard analysis the remaining compo- 
nents of the event rate computation are inputs taken 
from other sources. In particular, the fluxes of atmo- 
spheric neutrinos are taken from the results of numerical 
calculations, like those of Refs. |13ll2*] . which make use 
of a mixture of primary cosmic ray spectrum measure- 
ments, models for the hadronic interactions and simula- 
tion of particle propagation [xj. 

The attainable accuracy in the independent deter- 
mination of the relevant neutrino oscillation parameters 
from non-atmospheric neutrino experiments makes it pos- 
sible to attempt an inversion of the strategy in the atmo- 
spheric neutrino analysis: to use the oscillation param- 
eters as inputs in the data analysis in order to extract 
from data the atmospheric neutrino fluxes. 

There are several motivations for such direct deter- 
mination of the atmospheric neutrino flux from experi- 
mental data. First of all it would provide a cross-check 
of the standard flux calculations as well as of the size of 
the associated uncertainties (which being mostly theoret- 
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ical are difficult to quantify). Second, a precise knowl- 
edge of atmospheric neutrino flux is of paramount im- 
portance for high energy neutrino telescopes both 
because they are the main background and they are 
used for detector calibration. Finally, such program will 
quantitatively expand the physics potential of future at- 
mospheric neutrino experiments. Technically, however, 
this program is challenged by the absence of a generic 
parametrization of the energy and angular functional de- 
pendence of the fluxes which is valid in all the range of 
energies where there is available neutrino data. 

In this contribution wc present the first results on 
this alternative approach to the determination of the at- 
mospheric neutrino fluxes: we will determine these fluxes 
from experimental data on atmospheric neutrino event 
rates, using all available information on neutrino oscil- 
lation parameters, cross-sections, as well as information 
from the different sources of experimental and theoreti- 
cal uncertainty. The problem of the unknown functional 
form for the neutrino flux is bypassed by the use of 
neural networks as interpolants, since they allow us to 
parametrize the atmospheric neutrino flux without hav- 
ing to assume any functional behavior. Additional details 
on our approach can be found in 

Indeed the problem of the deconvolution of the at- 
mospheric flux from experimental data on event rates is 
rather close in spirit to the determination of parton dis- 
tribution functions in deep-inelastic scattering from ex- 
perimentally measured structure functions 0. For this 
reason, in this work we will apply for the determina- 
tion of the atmospheric neutrino fluxes a general strat- 
egy originally designed to extract parton distributions in 
an unbiased way with faithful estimation of the uncer- 
tainties 1 030. 



2 General strategy 

The general strategy that will be used to determine the 
atmospheric neutrino fluxes was first presented in Ref. [Hj 
(see also QHl)- It involves two distinct stages in order to 
go from the data to the flux parametrization. In the first 
stage, a Monte Carlo sample of replicas of the experi- 
mental data on neutrino event rates ( "artificial data" ) 
is generated. These can be viewed as a sampling of the 
probability measure on the space of physical observables 
at the discrete points where data exist. In the second 
stage one uses neural networks to interpolate between 
points where data exist. In the present case, this second 
stage in turn consists of two sub-steps: the determina- 
tion of the atmospheric event rates from the atmospheric 
flux in a fast and efficient way, and the comparison of 
the event rates thus computed to the data in order to 
tune the best-fit form of input neural flux distribution 

This strategy has also been successfully applied with dif- 
ferent motivations in oth er c ontexts like tau lepton decays 
[T7j and B meson physics ^Sj 



( "training of the neural network" ) . Combining these two 
steps, the space of physical observables is mapped onto 
the space of fluxes, so the experimental information on 
the former can be interpolated by neural networks in the 
latter. 

In the present analysis we use the latest data on at- 
mospheric neutrino event rates from the Super Kamiokande 
Collaboration ^ . Higher energy data from neutrino tele- 
scopes like Amanda |Hj is not publicly available in a for- 
mat which allows for its inclusion in the present analysis 
and its treatment is left for future work. 

The latest Super Kamiokande atmospheric neutrino 
data sample is divided in 9 different types of events: 
contained events in three energy ranges, Sub-GeV, Mid- 
GeV and Multi-GeV electron- and muon-like, partially 
contained muon-like events and upgoing stopping and 
througoing muon events. Each of the above types of 
events is divided in 10 bins in the final state lepton zenith 
angle . Therefore we have a total of Adat = 90 experi- 
mental data points, which we label as 

R^, i=l,...,JVdat. (1) 

Note that each type of atmospheric neutrino event rate 
is sensitive to a different region of the neutrino energy 
spectrum. In particular we note that for energies larger 
than E v ~ few TeV and smaller than E v ~ 0.1 GeV 
there is essentially no information on the atmospheric 
flux from the available event rates. 

The purpose of the artificial data generation is to pro- 
duce a Monte Carlo set of 'pseudo-data', i.e. A rep repli- 
cas of the original set of Ad a t data points, ^ art ^ fe ) such 
that the A rep sets of Adat points are distributed accord- 
ing to an Adat _ dimensional multi-gaussian distribution 
around the original points, with expectation values equal 
to the central experimental values, and error and covari- 
ance equal to the corresponding experimental quantities. 
This is achieved by defining 

Rt tm = Rt P) + rfM ot , k = 1, A rep , (2) 

where A rop is the number of generated replicas of the 

experimental data, and where are univariate gaus- 
sian random numbers with the same correlation matrix 
as experimental data. We can then generate arbitrarily 
many sets of pseudo-data, and choose the number of 
sets A rop in such a way that the properties of the Monte 
Carlo sample reproduce those of the original data set to 
arbitrary accuracy. 

In our case each neural network parametrizes a neu- 
trino flux. Since the precision of the available experimen- 
tal data is not enough to allow for a separate determi- 
nation of the energy, zenith angle and type dependence 
of the atmospheric flux, in this work we will assume the 
zenith and type dependence of the flux to be known with 
some precision and extract from the data only its energy 
dependence. That is, if NN (E v ) is the neural network 
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output when the input is the neutrino energy E v , then 
the neural flux parametrization will be 



^ (E v , c u ,t)= NN (E u ) ^ (E v , c u , t) 



(3) 



is a reference differential flux, which we take to be 
the most recent computations of either the Honda or 
the Bartol ^2. collaborations, which have been extended 
to cover also the high-energy region by consistent match- 
ing with the Volkova fluxes. 

Now let us describe a fast and efficient technique 
to evaluate the atmospheric neutrino event rates R* 1 
for an arbitrary input neural network atmospheric flux 
<pl nct ) [E v ,c v ,t). For example, for contained events any 
rate R,{ net ^ can be computed 



R 



(net) _ 



a,0,± 



oo r +l 



n t T dh dc » / dE v / dE; 



Emin J Bmin 



+ 1 



d,C n 



dE,, dc,, dh 



(E v ,c y ,h) 



P^ (E v ,c v ,h\r)) 



d 2 4 
dE t dc c 



{E v ,Ei,c a )e$ n [E u ci), (4) 



where (<£ Q ) is the flux of atmospheric neutrinos (an- 
tineutrinos) of type a, ctJ is the charged-current 

neutrino- (antineutrino-) nucleon interaction cross sec- 
tion, and see ^U] for a more detailed description of the 
above expression. The above expressions are very time- 
consuming to be evaluated, and this is specially a prob- 
lem in our case since in the neural network approach one 
requires a very large number of evaluations. 

However, the procedure can be speed up if one real- 
izes that for a given flux, the expected event rates can 
always be written as 

/l />oo 
dc v \ dE v ^{E v ,c v ,t)C i {E v ,c v ,t) , 
-1 

where Ci (E v ,c u ,t) contain all the E v and c v flux inde- 
pendent pieces in Eq.QJ. Consequently if one discretizes 
the differential fluxes in N e energy intervals and N z zenith 
angle intervals 



<z> ^ {E v , Cv ,t) = V v^fe (E v g J e ) o ( Cu g /,) 



e.z 



before the training, since they do not depend on the 
parametrization of the atmospheric neutrino flux. 

The determination of the parameters that define the 
neural network, its weights, is performed by maximum 
likelihood. This procedure, the so-called neural network 
training, proceeds by minimizing an error function, which 
coincides with the % 2 of the experimental points when 
compared to their theoretical determination obtained us- 
ing the given set of fluxes: 



({wj) = min € 



E^+E£ 



(8) 



i.thco 
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(W) 



R 



(art)(fe) 
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/ 



(9) 

The error function has to be minimized with respect to 
{wi}, the parameters of the neural network. The mini- 
mization of Eq. is performed with the use of genetic 
algorithms JS] . Because of the nonlinear dependence of 
the neural net on its parameters, and the nonlocal depen- 
dence of the measured quantities on the neural net (event 
rates are given by multi-dimensional convolutions of the 
initial flux distributions), a genetic algorithm turns out 
to be the most efficient minimization method. 

Thus at the end of the procedure, we end up with 
iV r ep fluxes, with each flux <£( not H fc ) given by a neural net. 
The set of N rcp fluxes provide our best representation 
of the corresponding probability density in the space of 
atmospheric neutrino fluxes: for example, the mean value 
of the flux at a given value of E v is found by averaging 
over the replicas, and the uncertainty on this value is the 
variance of the values given by the replicas. 



3 Results and discussion 

Now we discuss the results of our determination of the at- 
mosferic neutrino fluxes. Note that any functional of the 
neutrino flux (like for example event rates and its associ- 
ated uncertainty) can be computed taking the appropri- 
ate moments over the constructed probability measure 



(6) N rcp 

it is possible to write the theoretical predictions as a sum (jp (E u ))) = T |(£( nct )( fc ) ( E v 

of the elements of a bin-integrated flux table S'ilJ?^, rep k=i 



(10) 



p(nct) _ r<i ?7/( net ) 

"•% — 2^1 ^ezt^ezt J 
ezt 



(7) 



where the coefficients C* 2t , which are the most time- 
consuming ingredient, need only to be precomputed once 



In Fig. ^ we show our results for the atmosferic neutrino 
flux as compared with the Honda and Bartol compu- 
tations as well as with some direct measurements from 
Amanda |S] • In Fig. [21 we show that the effects of the 
uncertainties in the neutrino oscillation parameters are 
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Fig. 1 Comparison of the 1 - a uncertainty band of our Fig. 3 Comparison of the neural network results with exper- 
parametrization of the atmospheric neutrino flux with recent imental data for the Multi-GeV data sample, 
computations and AMANDA data. 
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Fig. 2 Results for fits varying the default neutrino oscilla- 
tion paramet ers from their central values within their 1-er 
uncertainties |14) . 



rather small with respect the intrinsic flux uncertain- 
ties. Finally, we show in Fig. [3] how our neural network 
parametrization successfully reproduces the features of 
experimental data. 

The work presented in this contribution can be ex- 
tended in several directions. First of all one could reduce 
the uncertainty in the large E v region by incorporating in 
the fit atmosferic neutrino measurements from neutrino 
detectors like Amanda [jJJ. Second, we could use at even 
larger energies simple parametrizations of the flux to ob- 
tain a better estimation of backgrounds at neutrino tele- 
scopes. Finally, we could estimate the required statistics 
that would be required in forthcoming atmosferic neu- 
trino experiments in order to determine from data also 
the zenith angle and flavor dependence of the neutrino 
flux. 
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